Wavefront depinning transition in discrete one-dimensional reaction-diffusion systems 
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Pinning and depinning of wavefronts are ubiquitous features of spatially discrete systems describ- 
ing a host of phenomena in physics, biology, etc. A large class of discrete systems is described by 
overdamped chains of nonlinear oscillators with nearest-neighbor coupling and controlled by con- 
stant external forces. A theory of the depinning transition for these systems, including scaling laws 
and asymptotics of wavefronts, is presented and confirmed by numerical calculations. 
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Spatially discrete systems describe physical reality in 
many different fields: propagation of nerve impulses 
along mielinated fibers pulse propagation through 

cardiac cells calcium release waves in living cells 
1^ , sliding of charge density waves Q , superconductor 
Josephson array junctions ||5|], motion of dislocations in 
crystals [||, atoms adsorbed on a periodic substrate [Q, 
arrays of coupled diode resonators [p| , and weakly cou- 
pled semiconductor superlattices A distinctive fea- 
ture of discrete systems (not shared by continuous ones) 
is the phenomenon of wavefront pinning: for values of a 
control parameter in a certain interval, wavefronts joining 
two different constant states fail to propagate Q. When 
the control parameter surpasses a threshold, the wave- 
front depins and starts moving The existence 
of such thresholds is thought to be an intrinsecally dis- 
crete fact, which is lost in continuum aproximations. The 
characterization of propagation failure and front depin- 
ning in discrete systems is thus an important problem, 
which is still poorly understood despite the numerous in- 
roads made in the literature [p]J^,^,|6|,pT 13 1. 

In this Letter, we study front depinning for infinite 
one-dimensional nonlinear spatially discrete reaction- 
diffusion (RD) systems. The nature of the depinning 
transition depends on the nonlinearity of the model, and 
is best understood as propagation failure of the travel- 
ing front. Usually, but not always, the wavefront profiles 
become less smooth as a parameter F (external field) de- 
creases. They become discontinuous at a critical value 
Fc- Below Fc, the front is pinned at discrete positions 
corresponding to a stable steady state. Fig. |l| shows 
wavefront profiles near the critical field. Individual points 
undergo abrupt jumps at particular times, which gives 
the misleading impression that the motion of the discrete 
fronts proceeds by successive jumps. Wavefront velocity 
scales with the field as \F — Fc\^ . For exceptional non- 
linearities, the wavefront does not lose continuity as the 
field decreases. In this case, there is a continuous tran- 
sition between wavefronts moving to the left for F > 
and moving to the right for F < Q: as for continuous sys- 
tems, front pinning occurs only at a single field F ~ Q. 



Wavefront velocity scales then linearly with the field. We 
discuss the characterization of the critical field (includ- 
ing analytical formulas in the strongly discrete limit), 
describe depinning anomalies (discrete systems having 
zero critical field), and give a precise characterization of 
stationary and moving fronts near depinning (including 
front velocity) by singular perturbation methods. Our 
approximations show excellent agreement with numeri- 
cal calculations. 

We consider chains of diffusively coupled overdamped 
oscillators in a potential V , subject to a constant external 
force F: 

diL 

^ Un+l-2Un+ Un-l+ F - Ag{un). (1) 

dt 

Here g{u) = V'{u) presents a 'cubic' nonlinearity, such 
that Ag{u) — F has three zeros, Ui < U2 < U3 in 
a certain force interval {g' {Ui{F / A)) > for i = 1,3, 
g'{U2{F/A)) < 0). Provided g{u) is odd with respect to 
112(0), there is a symmetric interval |F| < Fc where the 
wavefronts joining the stable zeros Ui{F/A) and U3{F/A) 
are pinned. For \F\ > Fc, there are smooth traveling 
wavefronts, u„(t) = u{n — ct), with u(— 00) — Ui and 
m(oo) = U3 The velocity c{A, F) depends on A 

and F and it satisfies cF < and c —> as |i^| —> i^c 
[|l5| . Examples are the overdamped Frenkel-Kontorova 
(FK) model (g = sinu) Mm and the quartic double well 
potential {V — {u^ — iy/4). Less symmetric nonlin- 
earities yield a non-symmetric pinning interval and our 
analysis applies to them with trivial modifications. 

Critical field. Stationary and traveling wavefronts can- 
not coexist for the same value of F [ [l5[. This fol- 
lows from a comparison principle for (p) |17]]. Pinning 
can be proved using stationary sub and supersolutions, 
which can be constructed provided the stationary so- 
lution is linearly stable. The largest eigenvalue of the 
linearization of (^ about a stationary profile Un{A,F), 
Un{t) = Un{A, F) + u„e^*, is given by 



\{A, F) = Min 



n+l 



v^f+Ag'{u,M.F))vl 



(2) 
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over a set of functions w„ which decay exponentially as 
n — > ±00. The critical field is uniquely characterized 
by \{A,Fc) = and A(A,F) < for \F\ < F^. Thus 
two facts distinguish the depinning transition: (i) one 
eigenvalue becomes zero, and (ii) stationary and moving 
wavefronts cannot coexist for the same values of the field. 

Equation shows that the critical field is positive 
for large A and typical nonlinearities. In fact, consider 
the FK potential. For F = there are two station- 
ary solutions which are symmetric with respect to U2, 
one taking on the value U2 (unstable dislocation), and 
the other one having m„ ^ U2 (stable dislocation) [^ . 
For large A, the stable dislocation has g'{un) > for 
all n, and (|) gives A(A,0) < 0. Since X{A,Fc) = 0, 
this implies that the critical field is nonzero. (A differ- 
ent proof can be obtained using the comparison principle 
[0,^). As A > decreases, several u„ may enter the 
region of negative slope g'{u): the number of points with 
g'{un) < increases as A decreases. It is then possible 
to have X{A, 0) = 0, i.e. F^ = 0, for a discrete system! 
Examples of this pinning anomaly will be given below. 

If F > 0, the stable dislocation is no longer symmetric 
with respect to U2- If F is not too large, all Un{A,F) 
avoid the region of negative slope g'{u) < 0. For large 
A and F and the generic potentials above mentioned, 
we have observed numerically that g' < for a single 
point, labelled uo{A,F). This property persists until Fc 
is reached. How does Fc depend on A7 For g = s'mu, 
it is well known that Fc vanishes exponentially fast as A 
goes to zero (the continuum limit). This was conjectured 
by Indenbom on the basis of a continuum approxi- 
mation, and numerically checked by Hobart JlSt in the 
context of the Peierls stress and energy for dislocations. 
More recently, Kladko et al derived the formula Fc — 
C exp(— TT^/ \J A — A'^/12) by means of a variational ar- 
gument. This argument can be used for other potentials 
and it suggests that Fc ~ Ce"''/^ as A ^ 0+ (with 
positive C and r/ independent of A) holds for a large class 
of nonlinearities King and Chapman have obtained 
an analogous result p3[ | using exponential asymptotics 

for the FK potential, Fc ^ C e^'^ 1 ^ and the wa vefront 
velocity after depinning, c ^ D ^ {F^ — F^)/A. This 
later result agrees with the scaling law c ^ \F — Fc\^ , 
found in a large class of discrete RD equations pT] , p^ . 

Anomalies of pinning. Despite widespread belief, it is 
not true that Fc > for all discrete systems. Using the 
characterization A(^, Fc) = 0, it is possible to see that 
having a zero critical field is equivalent to having a one- 
parameter family of continuous increasing stationary pro- 
files Un — u{n+a) satisfying u{x+l) + u{x — l) — 2u{x) — 
Ag{u{x)), with u(— 00) = J7i, u(oo) = C/3. In this case, a 
standard perturbation argument yields a wavefront speed 
having the same scaling as the continuum approxima- 
tion to the discrete system, c ~ C F as F 0. An 
example of nonlinearity presenting this anomalous pin- 
ning 1 20 1 can be obtained from u{x) = tanhx: it obeys 
the above equation with A = 1, Ui — —I, U3 = I and 



g{u) = -2tanh^(l)it(l-M2)/[l-tanh^(l)u2]. Further- 
more the wavefront velocity after depinning obeys the 
relation, c ~ — 3F/2 as F ^ 0. Thus nonlinearities pre- 
senting anomalous depinning belong to a different univer- 
sality class: the wavefront velocity has a critical exponent 
1 (and Fc ~ 0) instead of 1/2, which is the usual case for 
discrete RD systems (having Fc > 0). 

Asymptotic theory of wavefront depinning. We shall 
study the depinning transition in the strongly discrete 
limit v4 ^ 1, in which the structure of the wavefront is 
particularly simple. Firstly, consider the symmetric sta- 
tionary profile with Un ^ U2 for F = 0. The front profile 
consists of two tails with points very close to Ui and U^, 
plus two symmetric points uq, ui in the gap region be- 
tween Ui and U^. As F > increases, this profile changes 
slightly: the two tails are still very close to Ui{F/A) and 
U^{F/A). As for the two middle points, ui gets closer 
and closer to whereas wq moves away from Ui. This 
structure is preserved by the traveling fronts above the 
critical field: there is only one active point most of the 
time, which we can adopt as our uq. Then we can ap- 
proximate U-i ~ C/i, Ml ~ in (p, thereby obtaining 



duQ 



Ui 



U3 



-2uo- Ag{uo) + F. (3) 



This equation has three stationary solutions for F < Fc, 
two stable and one unstable, and only one stable sta- 
tionary solution for F > Fc. The critical field Fc is such 
that the expansion of the right hand side of (||) about the 
two coalescing stationary solutions has zero linear term, 
2 + Ag'{uo) = 0, and 



2uo + Ag{uQ) - Ui 



Us 



Fc 



(4) 



These equations for Fc and uo{A,Fc) have been solved 
for the FK potential, for which uq = cos~^(— 2/y4) and 
Ui + U;} = 2 sin^^(Fc/A)-|-27r. The results are depicted in 
Fig. and show excellent agreement with the numerical 
solution of (|^) for A > 10. Our approximation performs 
less well for smaller A, and it breaks down aX A — 2 with 
the prediction Fc = 0. Notice that Fc{A) ^ A as A in- 
creases. In practice, only steady solutions are observed 
for very large A. 

Let us now construct the profile of the traveling wave- 
fronts after depinning, for F sligthly above Fc. Then 
Uo(i) = 1*0(^7 Fc) + ^^o('^) obeys the following equation: 



^«(F-Fc)+A|5>o)||, 



(5) 



where we have used 2 -t- Ag'{uQ) — Q, ^ and ignored 
terms of order (F — Fc) /A and higher. This equation has 
the (outer) solution 



vo{t) 



l2{F-Fc) 
A\g"iuo)\ 



tan 



lA\g"iuo)\iF^Fc) 



{t - h) 

(6) 
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which is very small most of the time, but it blows up when 
the argument of the tangent function approaches ±7r/2. 
Thus the outer approximation holds over a time inter- 
val {t~to) - 71^/2/ y/A\g"{uo)\ {F - F^), which equals 
7r\/2(A2 - 4)-3(F - Fc)-i for the FK potential. The 
reciprocal of this time interval yields an approximation 
for the wavefront velocity, 

or |c| - {A^-A)i{F-F^)i/{'KV2) for a FK potential. In 
Fig. ^ we compare this approximation with the numeri- 
cally computed velocity for A ^ 100 and A ^ 10. 

When the solution begins to blow up, the outer so- 
lution ^ is no longer a good approximation, for uo(t) 
departs from the stationary value uo{A, Fc). We must go 
back to (|^) and obtain an inner approximation to this 
equation. As F is close to Fc and uo{t) — ua{A,Fc) 
is of order 1, we solve numerically (||) at F = Fc 
with the matching condition that uo{t) — uo{A, Fc) ^ 

2/[tt^\ A\g"{u„)\/{F - Fc) - A\g"{uo)\ {t - to% as [t - 
to) — > —oo. This inner solution describes the jump of uq 
to values close to U3. During this jump, the motion of 
Mo forces the other points to move. Thus, u-i{t) can be 
calculated by using the inner solution in (|l|) for uq, with 
F = Fc and u_2 « J7i. A composite expansion [£1| con- 
structed with these inner and outer solutions is compared 
to the numerical solution of (|l|) in Fig. ^. 

Notice that (j^) is the normal form associated with a 
saddle-node bifurcation in a one dimensional phase space. 
The wavefront depinning transition is a global bifurca- 
tion with generic features: each individual point u„(t) 
spends a long time, which scales as \F — Fd^i , near 
discrete values Un{A,Fc), and then jumps to the next 
discrete value on a time scale of order 1. The travel- 
ing wave ceases to exist for F < Fc- For these field 
values, discrete stationary profiles Un{A,F) are found. 
The above calculations give a normal form of the type 
(Pvo/dt^ = a{F — Fc) + f3vl instead of (||) for conserva- 
tive discrete systems (two time derivatives instead of one 
in (|l|)). The solution of this equation blows up in finite 
time as {F — Fc)^^, which gives a critical exponent of 
1/4 for the wavefront velocity near the critical field. 

The approximations to Fc{A) and the wavefront speed 
provided by the previous asymptotic theory break down 
for small A. In particular, for the FK potential and 
A < 2, no double zeroes of 2a; -|-^sin(a;) — (F-|-[/i-f?73) are 
found for F ^ Fc- What happens is that we need more 
than one point to approximate wavefront motion. Depin- 
ning is then described by a reduced system of more than 
one degree of freedom corresponding to active points. 
There is a saddle-node bifurcation in this reduced sys- 
tem whose normal form is of the same type as (^) . The 
jump of the active points after blow up is found by solving 
the reduced system with a matching condition p^ . As 
we approach the continuum limit, more and more points 



enter the reduced system of equations and exponential 
asymptotic methods become a viable alternative to our 
methods. 

In conclusion, we have studied depinning of wavefronts 
in discrete RD equations. The normal depinning tran- 
sition can be viewed as a loss of continuity of traveling 
front profiles as the critical field is approached: below the 
critical field, the fronts become pinned stationary profiles 
with discontinuous jumps at discrete values u„. In the 
strongly discrete limit, the critical field and these fronts 
can be approximated by singular perturbation methods 
which show excellent agreement with numerical solutions. 
The leading order approximation to the wavefront veloc- 
ity is then correctly given (scaling and prefactor) near 
the critical field. Depinning transitions for discrete RD 
equations apparently belong to two different universality 
classes. In the normal class, the wavefront velocity has 
a critical exponent 1/2. For certain nonlinearities, the 
stationary fronts are continuous functions of the discrete 
index at zero field. Then the critical field is zero, the de- 
pinning transition between stationary and moving fronts 
is continuous, with a critical exponent 1. This situation 
is the same as for continuous RD equations and we have 
called it anomalous pinning. 

AC thanks S. Hastings and J.B. McLeod for fruifful 
discussions. 
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FIG. 1. Traveling wavefront profiles near F — Fc for the 
FK potential and: (a) A = 2, (b) A = 100. 
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FIG. 2. (a) uo{A,Fc); (b) Critical field as a function of A; 
(c) Absolute and relative errors in Fc{A). 
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FIG. 3. Wavefront velocity versus F near F — Fc for the 
FK potential and: (a) A = 100, (b) A = 10. 
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FIG. 4. Wavefront profiles near Fc for the FK potential 
and A = 100. We show the three largest jumps in Fig. ^(b). 
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